ISOMETRIC IMMERSIONS AND COMPENSATED COMPACTNESS 



GUI-QIANG CHEN, MARSHALL SLEMROD, AND DEHUA WANG 

Abstract. A fundamental problem in diflterential geometry is to characterize intrinsic 
metrics on a two-dimensional Riemannian manifold which can be realized as iso- 
metric immersions into R^. This problem can be formulated as initial and/or boundary 
value problems for a system of nonlinear partial differential equations of mixed elliptic- 
hyperbolic type whose mathematical theory is largely incomplete. In this paper, we 
develop a general approach, which combines a fluid dynamic formulation of balance laws 
for the Gauss-Codazzi system with a compensated compactness framework, to deal with 
the initial and/or boundary value problems for isometric immersions in R"^. The com- 
pensated compactness framework formed here is a natural formulation to ensure the 
weak continuity of the Gauss-Codazzi system for approximate solutions, which yields 
the isometric realization of two-dimensional surfaces in R'^. 

As a first application of this approach, we study the isometric immersion problem 
for two-dimensional Riemannian manifolds with strictly negative Gauss curvature. We 
prove that there exists st C^'^ isometric immersion of the two-dimensional manifold in 

satisfying our prescribed initial conditions. To achieve this, we introduce a vanishing 
viscosity method depending on the features of initial value problems for isometric im- 
mersions and present a technique to make the apriori estimates including the L°° control 
and iif~^-compactness for the viscous approximate solutions. This yields the weak con- 
vergence of the vanishing viscosity approximate solutions and the weak continuity of the 
Gauss-Codazzi system for the approximate solutions, hence the existence of an isometric 
immersion of the manifold into R^ satisfying our initial conditions. 



1. Introduction 

A fundamental problem in differential geometry is to characterize intrinsic metrics on a 
two-dimensional Riemannian manifold which can be realized as isometric immersions 
into (cf. Yau [39]; also see [201 ESl |35]). Important results have been achieved for 
the embedding of surfaces with positive Gauss curvature which can be formulated as an 
elliptic boundary value problem (cf. [SO])- For the case of surfaces of negative Gauss 
curvature where the underlying partial differential equations are hyperbolic, the compli- 
mentary problem would be an initial or initial-boundary value problem. Hong in [22] first 
proved that complete negatively curved surfaces can be isometrically immersed in if the 
Gauss curvature decays at certain rate in the time-like direction. In fact, a crucial lemma 
in Hong [22] (also see Lemma 10.2.9 in [20]) shows that, for such a decay rate of the nega- 
tive Gauss curvature, there exists a unique global smooth, small solution forward in time 
for prescribed smooth, small initial data. Our main theorem, Theorem 5.1(i), indicates 
that in fact we can solve the corresponding problem for a class of large non-smooth initial 
data. Possible implication of our approach may be in existence theorems for equilibrium 
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configurations of a catenoidal shell as detailed in Vaziri-Mahedevan [38j . When the Gauss 
curvature changes sign, the immersion problem then becomes an initial-boundary value 
problem of mixed elliptic-hyperbolic type, which is still under investigation. 

The purpose of this paper is to introduce a general approach, which combines a fluid 
dynamic formulation of balance laws with a compensated compactness framework, to deal 
with the isometric immersion problem in (even when the Gauss curvature changes 
sign). In Section 2, we formulate the isometric immersion problem for two-dimensional 
Riemannian manifolds in via solvability of the Gauss-Codazzi system. In Section 3, we 
introduce a fluid dynamic formulation of balance laws for the Gauss-Codazzi system for 
isometric immersions. Then, in Section 4, we form a compensated compactness framework 
and present one of our main observations that this framework is a natural formulation to 
ensure the weak continuity of the Gauss-Codazzi system for approximate solutions, which 
yields the isometric realization of two-dimensional surfaces in M^. 

As a first application of this approach, in Section 5, we focus on the isometric immersion 
problem of two-dimensional Riemannian manifolds with strictly negative Gauss curvature. 
Since the local existence of smooth solutions follows from the standard hyperbolic theory, 
we are concerned here with the global existence of solutions of the initial value problem 
with large initial data. The metrics gij we study have special structures and forms usually 
associated with 

(i) the catenoid of revolution when gn = §22 = cosh(x) and gi2 = 0; 

(ii) the helicoid when gn = + y^, g22 = 1, and gi2 = 0. 

For these cases, while Hong's theorem [22j applies to obtain the existence of a solution 
for small smooth initial data, our result yields a large-data existence theorem for a C^'^ 
isometric immersion. 

To achieve this, we introduce a vanishing viscosity method depending on the features 
of the initial value problem for isometric immersions and present a technique to make the 
apriori estimates including the control and i?~^-compactness for the viscous approx- 
imate solutions. This yields the weak convergence of the vanishing viscosity approximate 
solutions and the weak continuity of the Gauss-Codazzi system for the approximate so- 
lutions, hence the existence of a C^'^-isometric immersion of the manifold into M'^ with 
prescribed initial conditions. 

We remark in passing that, for the fundamental ideas and early applications of compen- 
sated compactness, see the classical papers by Tartar [37] and Murat [H]. For applications 
to the theory of hyperbolic conservation laws, see for example [HIOl 112 ^ 117 ^ 1361. In particu- 
lar, the compensated compactness approach has been applied in [Sj EJ [ini ttH EH ES] to the 
one-dimensional Euler equations for unsteady isentropic flow, allowing for cavitation, in 
Morawetz \28\ [29] and Chen-Slemrod-Wang [7] for two-dimensional steady transonic flow 
away from stagnation points, and in Chen-Dafermos-Slemrod-Wang [5] for subsonic-sonic 
flows. 



2. The Isometric Immersion Problem for Two-Dimensional Riemannian 

Manifolds in 

In this section, we formulate the isometric immersion problem for two-dimensional Rie- 
mannian manifolds in via solvability of the Gauss-Codazzi system. 
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Let C be an open set. Consider a map r : ^ M'^ so that, for (x, y) G Q,, the two 
vectors {dxt, dyv} in span the tangent plane at r{x, y) of the surface r(Q) C M^. Then 

dxY X dyT 

n = — 

|a^r X dyr\ 

is the unit normal of the surface r(r2) C M'^. The metric on the surface in is 

ds^ = dr-dr (2.1) 

or, in local (rr, y)-coordinates, 

ds2 ^ (g^r . 5^ J.) + 2{dxY ■ dyT) dxdy + {dyV ■ dyr) {dyf. (2.2) 

Let gij,i,j = 1,2, be the given metric of a two-dimensional Riemannian manifold M 
parameterized on fi. The first fundamental form / for on Q is 

/ := gnidxf + 2gi2dxdy + g22{dyf- (2.3) 

Then the isometric immersion problem is to seek a map r : Q ^ M.^ such that 

dr ■ dr = /, 

that is, 

dxr-dxr = gu, d^r ■ dyV = g^, dyV ■ dyV = g22, (2.4) 
so that {dxY^dyv} in are linearly independent. 

The equations in (j2.4p are three nonlinear partial differential equations for the three 
components of r(x, y). 

The corresponding second fundamental form is 

// := -dn ■ dr = hii{dxf + 2hi2dxdy + h22{dyf , (2.5) 

and {hij)i<ij<2 is the orthogonality of n to the tangent plane. Since n • dr = 0, then 
d{n ■ dr) = implies 

-/I + n-d^r = 0, i.e., II = {n ■ d^r) (dx)"^ + 2{n ■ d^yr) dxdy + {n ■ d^r) {dy)"^ . 

The fundamental theorem of surface theory (cf. \13\ [20] ) indicates that there exists a 
surface in M'^ whose first and second fundamental forms are I and II if the coefficients 
(gij) and (hij) of the two given quadratic forms I and II with (gij) > satisfy the Gauss- 
Codazzi system. It is indicated in Mardare [27J (Theorem 9; also see [26]) that this theorem 
holds even when (/ijj) is only in for given {gij) in C^'^, for which the immersion surface 
is C^'^. This shows that, for the realization of a two-dimensional Riemannian manifold 
in with given metric {gij) > 0, it suffices to solve (hij) S L°° determined by the 
Gauss-Codazzi system to recover r a posteriori. 

The simplest way to write the Gauss-Codazzi system (cf. [HI [20]) is as 



dxM - dyL = vf^L - 2Tf^M + vf^N, 
dxN - dyM = -r^2^L + 2V^l^M - T^^^N, 



(2.6) 



with 

LN-M"^ = K. (2.7) 

Here 

h^^ _ _ h^o _ _ /199 
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\g\ = det{gij) = 511522 — 512) i^i^iU) is the Gauss curvature that is determined by the 
relation: 

^(^ n,\ — ^^^^^ p — r, f f> r("^) f> r(™) _i_ r(")r(™) r(")r(™)^ 
~ "iffT' ~ ^ ~ + ^ ^ nfc - ifc ^ nj j ' 

Rijki is the curvature tensor and depends on (gij) and its first and second derivatives, and 

^if = {dm + diQji - digij) 

is the Christoffel symbol and depends on the first derivatives of {gij), where the summation 
convention is used, (5^') denotes the inverse of (gij), and (81,82) = {8x,8y). 

Therefore, given a positive definite metric {gij) € C^'^, the Gauss-Codazzi system gives 
us three equations for the three unknowns (L, Af , N) determining the second fundamental 
form //. Note that, although {gij) is positive definite, -R1212 may change sign and so 
does the Gauss curvature k. Thus, as we will discuss in Section 3, the Gauss-Codazzi 
system (I2.6p -( f2?7|l generically is of mixed hyperbolic-elliptic type, as in transonic flow (cf. 
[21 13 El EQ])- In §3-4, we introduce a general approach to deal with the isometric immersion 
problem involving nonlinear partial differential equations of mixed hyperbolic-elliptic type 
by combining a fluid dynamic formulation of balance laws in §3 with a compensated 
compactness framework in §4. As an example of direct applications of this approach, in 
§5, we show how this approach can be applied to establish an isometric immersion of a 
two-dimensional Riemannian manifold with negative Gauss curvature in M^. 

3. Fluid Dynamic Formulation for the Gauss-Codazzi System 



From the viewpoint of geometry, the constraint condition (j2.7p is a Monge-Ampere 
equation and the equations in (|2.6p are integr ability relations. However, our goal here 
is to put the problem into a fluid dynamic formulation so that the isometric immersion 
problem may be solved via the approaches that have shown to be useful in fluid dynamics 
for solving nonlinear systems of balance laws. To achieve this, we formulate the isometric 
immersion problem via solvability of the Gauss-Codazzi system (12. 6p under constraint 
()2.7p . that is, solving first for hij,i,j = 1,2, via ()2.6p with constraint ()2.7p and then 
recovering r a posteriori. 

To do this, we set 

L = pv^ +p, M = —puv, N = pu^ + p, 

and set = v? + as usual. Then the equations in (|2.6p become the familiar balance 
laws of momentum: 

8x {puv) + dy {pv" + p) = - {pv" + p)rg - IpuvV^^l - {pu^ + p)Tf} , ^^^^ 

8x{pu^ +p)+ dy{puv) = -{pv^ +p)r^y - 2puvV^ll - (pn2 +p)rS\\ 

and the Monge-Ampere constraint (12. 7p becomes 

PV(^ + = K. (3-2) 

From this, we can see that, if the Gauss curvature k is allowed to be both positive and 
negative, the "pressure" p cannot be restricted to be positive. Our simple choice for p is 
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the Chaplygin-type gas: 

1 

P = — • 
P 



2 1 

-q + ^ = 



Then, from ()3.2p . we find 



and hence we have the "Bernouhi" relation: 

P=^4=- (3-3) 

This yields 

P = -\/q^, (3.4) 

and the formulas for u'^ and v^: 

=p{p- M), v'^=p{p-L), M"^ = {N -p){L -p). 

The last relation for gives the relation for p in terms of (L, M, N)^ and then the first 
two give the relations for {u,v) in terms of {L, M, N). 
We rewrite (13.11) as 



dxipuv) + dy{pv'^ +p) = Ri, 
dx{pv? +p) + dy{puv) = R2, 



2 , ^ . . . ^ ^ (3-5) 



where Ri and R2 denote the right-hand sides of (j3.ip . 

We now find the corresponding "geometric rotationality-continuity equations" . Multi- 
plying the first equation of ()3.5p by v and the second by u, and setting 

dxV — dyU = —a, 



we see 



—div( pu, pv) dyK = — - + au, 

P 2 p 

-div{pu,pv) - l-dxK. = — - av, 
P 2 p 



and hence 



(3.6) 



dw{pu,pv) = --dyK-\ \ , 

dw{pu,pv) = --dxK-\ . 

2u u u 

Thus, the right hand sides of (j3.6p are equal, which gives a formula for a: 

0- = -^{^C^PC>xi^ + R2) - ui^pdyK + Ri)^ . (3.7) 
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dxV - dyU = -z^[u{-pdyK + Ri) - v{-pdxK + ^^2) ), (3.8) 
dxipu) + dy{pv) = ^^d^K + ^^dyK + -^Ri + -^i?2- (3.9) 



If we substitute this formula for a into (j3.6p . we can write down our "rotationality- 
continuity equations" as 

l^(u{^pdyK + Ri) -v{^, 

1 pu 1 o ^ ^ 

-^dxK + -^OyK + -^Ri + — 

2 2 q^ ^ q^ q^ 

In summary, the Gauss-Codazzi system (|2.6p - ()2.7p . the momentum equations (|3.ip - (j3.4p . 
and the rotationahty-continuity equations (j3.3p and (|3.8p - (j3.9p are all formally equivalent. 

However, for weak solutions, we know from our experience with gas dynamics that this 
equivalence breaks down. In Chen-Dafermos-Slemrod-Wang [5], the decision was made 
(as is standard in gas dynamics) to solve the rotationality-continuity equations and view 
the momentum equations as "entropy" equalities which may become inequalities for weak 
solutions. In geometry, this situation is just the reverse. It is the Gauss-Codazzi system 
that must be solved exactly and hence the rotationality-continuity equations will become 
"entropy" inequalities for weak solutions. 

The above issue becomes apparent when we set up "viscous" regularization that pre- 
serves the "divergence" form of the equations, which will be introduced in §5.3. This is 
crucial since we need to solve (|3.8p - (|3.9p exactly, as we have noted. 

To continue further our analogy, let us define the "sound" speed: 

c^=p'{p), (3.10) 

which in our case gives 

= 1. (3.11) 



P 

Since our "Bernoulli" relation is (13.31). we see 



= q^ + K. (3.12) 

Hence, under this formulation, 

(i) when k > 0, the "flow" is subsonic, i.e., q < c, and system (j3.ip ~ (|3.2p is elliptic; 

(ii) when k < 0, the "flow" is supersonic, i.e., q > c, and system ()3.ip - (j3.2p is hyperbolic; 

(iii) when k = 0, the "flow" is sonic, i.e., q = c, and system ()3.ip - (|3.2p is degenerate. 

In general, system ()3.ip - ()3.2p is of mixed hyperbolic-elliptic type. Thus, the isomet- 
ric immersion problem involves the existence of solutions to nonlinear partial differential 
equations of mixed hyperbolic-elliptic type. 

4. Compensated Compactness Framework for Isometric Immersions 

In this section, we form a compensated compactness framework and present our new 
observation that this framework is a natural formulation to ensure the weak continuity of 
the Gauss-Codazzi system for approximate solutions, which yields the isometric realization 
of two-dimensional Riemannian manifolds in M^. 

Let a sequence of functions (L^, M^, A^^)(x, y), defined on an open subset C M^, 
satisfy the following Framework (A): 

(A.l) |(L^,M^, A^^)(x,y)| < C a.e. (x,y) G 0,, for some C > independent of e; 
(A. 2) dxM^ — dyU and dxN^ — dyM^ are confined in a compact set in ff;~^(i7); 
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(A. 3) There exist o^(l), j = 1, 2, 3, with o^(l) ^ in the sense of distributions as e ^ 
such that 

d^M' - dyL' = T^^L' - 2rgM^ + vf-^ + o^(l), 

(4 1) 

d^N' - dyM' = -r^L^ + ir'l^M' - T^^^N^ + o|(l), 

and 

L^N^ -{M^f = K + ol{l). (4.2) 

Then we have 

Theorem 4.1 (Compensated compactness framework). Let a sequence of functions 
{L^ ,N^){x,y) satisfy Framework (A). Then there exists a subsequence (still labeled) 
{L^ , , N^){x,y) that converges weak-star in L°°{Q) to {L,M,N) as e — > such that 

(i) \{L,M,N){x,y)\ <C a.e. (x,y) gJ]; 

(ii) the Monge-Ampere constraint (|2.7p is weakly continuous with respect to the sub- 
sequence {L^ , M"^ , N^)(x,y) that converges weak-star in L°°{i}) to {L,M,N) as 

(iii) the Gauss-Codazzi equations in (|2.6p hold. 

That is, the limit (L, M, N) is a bounded weak solution to the Gauss-Codazzi system (j2.6p - 
()2.7p . which yields an isometric realization of the corresponding two-dimensional Riemann- 
ian manifold in R'^ . 



Proof. By the div-curl lemma of Tartar-Murat [371 [3T] and the Young measure repre- 
sentation theorem for a uniformly bounded sequence of functions (cf. Tartar [37]), we 
employ (A.1)-(A.2) to conclude that there exist a family of Young measures {vx,y}[x,y)&u 
and a subsequence (still labeled) {L^ , , N^){x,y) that converges weak-star in L°°{Q) to 
(L, M,N) as e ^ such that 

(a) {L,_M,N){x,y) = {{ux,y, L),{vx,y, M),{ux,y, N)) a.e. {x,y) G Q; 

(b) \{L,M,N){x,y)\<C a.e. (x,y) Gfi; 

(c) the following commutation identity holds: 

{i^x,y, - LN) = {vx,y, Mf - {ux,y, L) {ux,y, N) = {Mf - LN. (4.3) 

Since the equations in (j4.ip are linear in [L^ ,N^), then the limit {L,M,N) also 
satisfies the equations in (j2.6p in the sense of distributions. 
Furthermore, condition (j4.2p yields that 

{ux,y, LN - M^) = k{x, y) a.e. (x, y) G (4.4) 

The combination (j4.3p with (j4.4p yields the weak continuity of the Monge-Ampere con- 
straint with respect to the sequence {L^,M^,N^) that converges weak-star in L°°(il) to 
{L,M,N) as e ^ 0: 

LN - {Mf = K. 

Therefore, {L, M, N) is a bounded weak solution of the Gauss-Codazzi system (I2.6p - (l2.7p . 
Then the fundamental theorem of surface theory implies an isometric realization of the 
corresponding two-dimensional Riemannian manifold in . This completes the proof. □ 
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Remark 4.1. In the compensated compactness framework, Condition (A.l) can be relaxed 
to the following condition: 

(A.l)' ||(L^M^ Af^)lliP(Q) <C, p>2, for some C > independent of e. 

Then all the arguments for Theorem 4.1 follow only with the weak convergence in LP(i}),p > 
2, replacing the weak-star convergence in L°°(J7), with the aid of the Young measure rep- 
resentation theorem for a uniformly bounded sequence of functions (cf. Ball ^). 

There are various ways to construct approximate solutions by either analytical methods, 
such as vanishing viscosity methods and relaxation methods, or numerical methods, such 
as finite difference schemes and finite element methods. Even though the solution to 
the Gauss-Codazzi system may eventually turn out to be more regular, especially in the 
region of positive Gauss curvature k > 0, the point of considering weak solutions here is to 
demonstrate that such solutions may be constructed by merely using very crude estimates. 
Such estimates are available in a variety of approximating methods through basic energy- 
type estimates, besides the -L°° estimate. On the other hand, in the region of negative 
Gauss curvature «; < 0, discontinuous solutions are expected so that the estimates can be 
improved at most up to BV in general. 

The compensated compactness framework (Theorem 4.1) indicates that, in order to 
find an isometric immersion, it suffices to construct a sequence of approximate solutions 
(L'^ , , N^){x,y) satisfying Framework (A), which yields its weak limit {L,M,N) to be 
an isometric immersion. To achieve this through the fiuid dynamic formulation ()3.ip and 
(|3.3p (or (|3.4p ). it requires a uniform estimate of {u^,v^) such that the sequence 

with 

1 



= = V(U=)2 + {v^y + K 



satisfies Framework (A). 

The fluid dynamic formulation, (j3.ip and (13. 3p (or (|3.4p ). and the compensated compact- 
ness framework (Theorem 4.1) provide a unified approach to deal with various isometric 
immersion problems even for the case when the Gauss curvature changes sign, that is, for 
the equations of mixed elliptic-hyperbolic type. 

5. Isometric Immersions of Two-Dimensional Riemannian Manifolds with 

Negative Gauss Curvature 

As a first example, in this section, we show how this approach can be applied to establish 
an isometric immersion of a two-dimensional Riemannian manifold with negative Gauss 
curvature in R^. 

5.1. Reformulation. In this case, k < in and, more specifically, 

K = —7^, 7 > in J7. 
For convenience, we assume 7 G C"^ in this section and rescale (L, M, N) in this case as 

L M ~ N 

L = -, M = — , = — , 

7 7 7 
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SO that (12.71) becomes 



LN -M'^ = -1. 

Then, without ambiguity, we redefine the "fluid variables" via 

L = pv^ + p, M = —puv, N = pu^ + p, 

and set = + v"^ , where we have still used {u,v,p,p) as the scaled variables and will 
use them hereafter (although they are different from those in §2-§4). 

Then the equations in (j2.6p become the same form of balance laws of momentum: 



dx{puv) + dy{pv^ +p) = Ri, 



,2 



dr^{pu + + dy{puv) = R2, 



(5.1) 



where 



i?i := -(pt;2 + p)fg) - 2puvff^ - {pu^ + p)tf} , (5.2) 

R2 := -{Pv^ + P)f g - 2puvt^S - {P^^ + P)f , (5.3) 

-•^ 11 — -"^ 11 "T ^ ' -"^ 12 — ^ 12 "T 2^ ' 22 — 22 ' 

^(2) _ -p(2) ^(2) _ ^(2) 7x -p{2) _ p(2) 7j/ 

11 — 11 ' 12 — 12 "T 2^ ' 22 — -"^ 22 "T ^ • 

Furthermore, the constraint LN — = — 1 becomes 

ppq^+p^ = -l. (5.4) 
From p = and (15. 4p . we have the "Bernoulli" relation: 

or p = —\/ — 1^ (5.5) 



which yields 

u^=p{p-N), v'^=p{p-L), {Mf = {N -p){L -p). (5.6) 

Then the last relation in (j5.6|) gives the relation for p in terms of {L,M, N), and the first 
two give the relations for (u, v) in terms of (L, M, A^). 

Similarly to the calculation in §3, we can write down our "rotationality-continuity 
equations" as 

dxV — dyU = ^ (vR2 — uRi) =: Si, (5-7) 

d^pu) + dyipv) = ^Ri + 4^2 =: ^2. (5.8) 
Under the new scaling, the "sound" speed is 

1 

P 

Then the "Bernoulli" relation (j3.3p yields 

= q^ - 1. (5.10) 
Therefore, q > c, and the "fiow" is always supersonic, i.e., the system is purely hyperbolic. 



2 - p'{p) = ^ > 0. (5.9) 
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5.2. Riemann invariants. In polar coordinates {u,v) = {qcos9,qsm9), we have 
Ri = pq^ cos^ 9 rfj - 2pq^ sin 9 cos 9 tf^^ + pq^ sin^ 9 ff^^ - p(f + ff^^) , 

i?2 = pq"^ cos^ 9 r^22 - 2pg^ sin 9 cos 9 t^^j + pq^ sin^ 9 fj^^ - p{V^^^ + fj\^) , 
and then (j5.7p and (jS.Sp become 

sin^5a;g + qcos9 dx9 — cos9dyq + (/sin^ = Si, (5-11) 



cos 6^ . a a n , ^mt^ 

—^5 rrCaig + sm9 dx9 + TrOyq - cos 6* dj^e' = 

That is, as a first-order system, ()5.7p and (|5.8p can be written as 



S2. (5.12) 



sin( 

1 

g(g^-i) 



(7 cos ( 
cos sin ^ 



9x 



+ 



cos 



qsmt 



7— TT sin 9 — cos 



a, 



'52 



. (5.13) 



One of our main observations is that, under this reformation, the two coefficient matrices 
in ()5.13p actually commute, which guarantees that they have common eigenvectors. The 
eigenvalues of the first and second matrices are 



X± = sin 9 ± 



cos ( 



cos 9 ± 



sm( 



and the common left eigenvectors of the two coefficient matrices are 

1 



(±- 



1 . 



qy^q'^ - I 

Thus, we may define the Riemann invariants W± = W±{9, q) as 

1 



deW± = 1, dgW± = ± 



which yields 



,1. 



W± = 9 ± arccos (-) . 



(5.14) 



(5.15) 



Now multiplication (|5.13p by (dqW±,deW±) from the left yields 



/ 2 _ 1 

A+ {dqW+ dxq + deW+ dj) + ^Ji+ {dqW+ dyq + deW+ dy9) = SidqW+ - ^ ^2, 



/ 2 _ 1 

A_ {dqW^ d^q + deW. dj) + /i_ {dqW- dyq + deW^ dy9) = SidqW^ - ^ 52 



Prom (I5T3D . 

dxW± = dgW+ dxq + dew+ dj, dyW± = dgW+ dyq + deW+ dy9, 



(5.16) 

52. 
(5.17) 
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then we can write (j5.16p and (j5.17|) as 
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\+d,W+ + ^i+dyW^ 



1 



1 



Jl '-'2, 



1 



(5.18) 
(5.19) 



5.3. Vanishing viscosity method via paraboHc regularization. Now we introduce a 
vanishing viscosity method via parabohc regularization to obtain the uniform L°° estimate 
by identifying invariant regions for the approximate solutions. 

First, if Ri and R2 denote the additional terms that should be added to the right-hand 
side of the Gauss-Codazzi system (15. ip . our first choice is 

(5.20) 



Ri = edyipv), R2 = edyipu), 

which gives us the system of "viscous" parabolic regularization: 

d^ipuv) + dyipv^ +p) = Ri+ edyipv) = Ri + R^, 
dx{pv? +p) + dy{puv) = R2 + £dy{pu) = R2 + R2- 

From equations (|3.8p and (jS.Op . we see 

= ~~2 (^-^2 - uRi^ , §2 = [vRi + uR 

should be added to Si and S2 on the right-hand side of (|3.8p and (|3.9p . In polar coordinates 
{u,v) = {q cos 6 , q sin 6) , (|5.22p becomes 

Si = e—dyOdyipq) + edle, §2 = e-dl{pq) - epidyOf 
pq q 



(5.21) 



(5.22) 



(5.23) 



Note the identity 



edy(^&Tccos (~)) = (arccsc(/3g)) 



-ed,. 



1 



pqV p^q^ - 1 
1 



pqVp^q 



2r,2 



1 



dyipq) 



P 
S2 
P' 



dyipq) - -^dl{pq) 



-{dyOf. 



Then 



and thus 



S2 

p2 



-edl 



arccos 



pqyjfp^ 



^)a,(pg)-i(a, 



2 



Si - (g2 - 1)^2 = Si-^ 

P^ 

= ~^y^ dyiPl) + + ^9y ( arccos (-)^ -|- edy 



PQVf^ 



dyipq) + -idy 
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Since 
then 

Si - {q^ - 1)52 = edlW+ + '^dyW+ dyipq) + e{dyW+)'' 

Similarly, using 

dyipq) 



dyO = dyW- 



pq\J p^tp' - l' 
we have 

-Si - (g2 - 1)^2 = -edlW. - "^dyW. dyipq) + eidyW.)\ 



Thus, if we add the above 5i and 5*2 to the original Si and 52, (|5.18p and (|5.19p become 
/ dW+ dW^^ 



dy 

edlW+ + ^dyW+ dyipq) + eidyW+)'' + 5i - {q^ - l) 52 



(5.24) 



dy 



(5.25) 
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(5.26) 



= -ed^yW^ - '^dyW- dyipq) + eidyW^f -Si- iq^ - l) 52 

Plugging Ri and i?2 into 5i and 52 yields 
Si ± iq^ - 1) 52 

= -q sin e (r'^j cos^ 6 - IT^j sin cos + tf^ sin^ - ^ (f^2^ + f 

- fg^ cos^ + 2ffj sin cos - ff^^ sin^ 6 + ^ (fg^ ' i^^^^ 
± i I g cos 6* (^f^a^ cos^ - 2f sin 9 cos 9 + f sin^ 9-^ (f + fj^^ 

+g sin 9 (fg^ cos^ - 2f [J sin cos + f sin^ 9-^ (fg^ + f 

Then system ()5.24p - ()5.25p is parabolic when A+ > and A_ < 0. 

Furthermore, setting iE,F,G) = ign, gi2, 922), we recall the following classical identi- 
ties: 

p(l) _ GE^-2FF^+FEy ^(1) _ 2GF^ - GG^ - FGx 

Ml- 2(EG-F-^) ' ^22- 2(SG-F2) 

p(2) _ 2EF^-EEy-FE^ p{2) _ EGy - 2FFy + FGx 

Ml- 2(EG-F^) ' ^22- 2iEG-F^) 

_ GEy-FG^ p{2) _ EGx - FEy 

12 - 2{EG-F-^) ' 12 - 2(^(^_^2)' 
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{EG - F^fn = det 



and 7^ = — k. 



1 



xy 

ir 



- K F 
F 



F 
G 



2^y 



det 





1 P 

2^y 

In 

,2^x 



i p 

2^y 
E 

F 



\Gx 

F 
G 



5.4. L°°— estimate for the viscous approximate solutions. Based on the calculation 
above for the Riemann invariants, we now introduce an approach to make the estimate. 
First we need to sketch the graphs of the level sets of W± . If 

W± = lb arccos (-) = C± for constants C±, 
9 



then 



dQ d ( A 
— = =F ~r arccos I - 
dq dq\ q 



q^yq'^ -1 



and, as g — > cxD, 
dO 



0, 



C± =F arccos (0) = C± =F 



vr 



dq - ' ' 2 

See Fig. [U for the graphs of the level sets W± = C± . 



on W± = G±, 



on W± = C±. 






V 




















u 




g= l\ 





Figure 1. Level sets 



Next we examine the meaning of inequality Wj^ < (7+, i.e. 

+ arccos (-) < C+. 



For example, \i q = 1, then 6 < C_|_. This indicates the region of W+ < C_|_ as sketched in 
Fig. [TJa). Similarly, W- > C_ means 

9 — arccos (-) > C_, 

q 
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and, if g = 1, then 9 > C_, and the region of > C_ is sketched in Fig. mb). Thus we 
see that 

W+ < C+ means the region below W+ = C+, 

W- > C_ means the region below W- = C^; 

and 

W+ > C7+ means the region above W+ = C+, 

W- < C- means the region above W- = C_. 

As an example, we now focus on the case that 

F = 0, E{x) = G{x). (5.27) 



Then 



Therefore, we have 

f(l) _ ^' , 7' -p(l) _ _ . _ -p(2) _ E' i ~(2) _ 

^11 -2^+^' ^12 -U, i22 - ^11-^' ^12-2£; + 27' ^22 -U, 

and the right-hand side of ()5.26p is equal to 

1 / k' 2 ^ 1 / k' 2^'\ 

^ ( A " ^ '^^^ J ^ ^ V 7 + ^ J 

Thus, the two solutions 9±{q) that make the right-hand side of ()5.26p equal to zero satisfy 

tan 6 = — (5.28) 

If we fix the intersection point of 0±{q) at 

61 = 0, q = qQ = p, 

where /? > 1 is a constant, then the above ordinary differential equation (j5.28p becomes 

1 k'(x) E'ix) 

+ "FTT = 0' (5-29) 



I.e., 

Thus, 

Since nix) < 0, then 



/32 k(x) ^(x) 

^£;(x)) 



-^ln(k(x)|^ 
ax 



|k(x)|^£'(x) = const. > 0. 



k(x) = —kqE{x) ^ , 
where kq > is a constant, and equation (j5.28p for tan^ becomes 



/32 - (/?2 _ i),2 • 

Assume that we have a solution £'(x) to (j5.29p . Fix another constant 1 < a < (3. Then 
the curves 9±{q) are independent of (x,y) and look like the sketch in Fig. [2j 
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V 




















j 


\ 

a 


/ /A 


u 













Figure 2. Graphs of Q± 

W- = W-{0,a} 

























J 


a \ 




u 






8-{q) 







,W+ = W+{0,a) 

Figure 3. Invariant regions 

Next, note that W±, when evaluated at 9 = 0,q = a, and 9 = 0,q = P, take on the 
constant values, i.e., independent of Equations (|5.24p and (|5.25p imply that, at any 

point where VW+ = (respectively VW- = 0), 



. j>0 foT9>9+{q), 



ejd;W. 



<Q ioT9<9+{q), 
>0 ioi9>9-{q), 



[<0, ioi9<9-{q). 
Hence, when > and A_ < 0, by the maximum principle (cf. [1 

has no internal maximum for 9 > 9j^{q), 
W+ has no internal minimum for 9 < 9^{q), 
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W- has no internal maximum for > 9-{q), 
W- has no internal minimum for 9 < 9-{q). 

Define 

W±(0,/3) = icos^i (i), W±(0,a) = icos-^ (^). 

Therefore, if the data is such that W+ < VF+(0, then Wj^ can have no internal maximum 
greater than VF+(0, /?) for 9 > 9^{q). Similarly, if the data is such that W- > W-{0,l3) 
(= VF+(0,/?)), then W- can have no internal minimum less than W-{0,f3) for 9 < 9-{q). 
Furthermore, if the data is such that W+ > VF4.(0,a), then W+ can have no internal 
minimum less than W+(0, a) for 9 < 9+{q); if that data is such that < W-{0,a), 
then W- can have no internal maximum greater than W-{0,a) for 6 > 9^{q). Thus, the 
diamond-shaped region in Fig. [3] provides the upper and lower bounds for W±. 
From the definition of X±: 

cos 6* 

A± = smt' ± — , 
we easily see that the lines A± = sketched in Fig. [H 




Figure 4. Graphs of A± = 

Notice that A+ > and A_ < are in the region above A+ = and below A_ = 0. If we 
now super-impose Fig. [3] on top of Fig. |4]and choose a sufficiently close to (3, we see that 
there is a region where the four-sided region of Fig. [3] is entirely confined in the region 
above A+ = and below A_ = in Fig. [H Hence, the parabolic maximum/minimum 
principles apply and the four-sided region is an invariant region. For example, in the 
half-plane: 

n := {{x,y) : x>0, y S M} (5.30) 
with periodic initial data {q{0,y), 9(0, y)) prescribed in the four-sided region, the maxi- 
mum/minimum principles yield the invariant region for the periodic solution. 

There is an alternative symmetry about ^ = f • If we set = V' + f j then the right-hand 
side of ()5.26p becomes 
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If i/;^ satisfy the above equations, then ip^ and are symmetric about 1/^ = 0, i.e., 6 = ^. 
Look for the crossing on ■(/; = so that 

1 / k' 2 E' 



With the crossing at g = /?, this gives 

-Ik' E' 



K E 

that is, 



+ — = 0, (5.31) 



\k{x)\~^ E{x) = const. 
To exploit the symmetry, we now take x as a space-hke variable and y as a time-like variable 
and replace dy by d'^. Now it is the interior between the lines = and /i+ = that gives 
the preferred signs /i+ > and /i_ < 0, which keep (I5.24p and (I5.25[) parabolic. Hence, 
similar to (i), in the case that the initial data and the metric E{x) = G{x), F{x) = 0, 
is periodic in x, then the periodic solution will stay in the four-sided invariant region in 
which the initial data lies. 

All these arguments yield the uniform L°° bounds for (n"^, t;'=,p^, p^), which implies 

\{L',M',N')\ < C, 
for some constant C > depending only on the data and ||7||ioo. 

Finally, let us examine the following examples: 

2 



Example 5.1. Catenoid: E{x) = {cosh{cx))^^^ , k{x) = —kqE{x)~^ , where c 7^ and 
Ko > are two constants. Substitution them into (|5.29p (where x is time-like and y is 
space-like) yields 

P>1. 

Of course, (j5.29p is satisfied when 



E{x) = {cosh{cx)f^^ k{x) = -kqE{x) 

with (3 > 1 (where x is space-like and y is time-like). 
Example 5.2. Helicoid: The metric associated with the helicoid is 

ds'^ = E{dXf + {dYf 
with E{Y) = + and the Gauss curvature 

where X > is a constant. To apply our previous result, for the special case of isothermal 
coordinates given in ()5.27p . we first make a change of variables to rewrite the helicoid 
metric in isothermal coordinates, i.e., allow X,Y to depend on x,y so that 

ds"^ = {EXl + Y^)dx'^ + 2{EX^Xy + Y^Yy)dxdy + {EX^ + Y^)dy\ 

Hence, if we set 

Y, = -VEXy, Yy = Vex,, 
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then 

ds^ = E{Xl + Xl){dx^ + dy"^), 

which gives the metric in isothermal coordinates. The above equations for X and Y may 
be rewritten as 

Y — Y 



/E " VE 

and with 



^=_ = ln(y+yA^), 



we have 



i.e., the Cauchy-Riemann equations. A convenient solution is given by 

4>=-x, X = y, 

which yields 

-x = in(y + VA2 + y2), 

that is, 

y = -i(AV-e--). 

Thus, in the new {x,y)- coordinates, 
and 

ds^ = (iA2 + 1(aV^ + {dx^ + dy^). 

Hence, we have 

2E'{x) _ k'{x) 
E{x) k{x) 
and so relation (|5.29|) is satisfied with (5 = \f2. 

Example 5.3. Torus: The metric for the torus is usually written as 

ds^ = EdX^ +b'^dY'^, 

with 

E=ia + bcosYf, <Y)= 

o(a + ocos Y) 

where a > 6 > are constants. The same argument as given in Example 5.2 above yields 
the metric in isothermal coordinates as 

ds"^ = E{dx^ + dy^) 

with 

E = {a + bcosYf = (a + 6cos((^"^(x))2), 

where 

b I \J c? — b'^ sin Y \ 



- 62 I b + acosY 
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, , cosy cosidi ^(x)) 

K[X) — — 



b{a + bcosY) b{a + bcos{<f>-^{x))) ' 
A direct computation yields 

k'{x) _ a(0-i(x))'tan(0-i(x)) E'{x) _ 6(0-i(x))' sin(0-i(x)) 
k{x) a + b cos{4>~^ (x)) ' E{x) a + bcos{(j)^^{x)) ' 

and the ratio 

k'{x) lE'ix) _ a 

k{x) I E{x) 6cos((/)~^(x)) 
is not a constant. So (|5.29p does not hold. Hence our Proposition 5.1 will not directly 
apply to that piece of the torus possessing negative Gauss curvature. 

5.5. if;^^— compactness. We now show how the //^^^-compactness can be achieved for 
the viscous periodic approximate solutions via parabohc regularization. 

In §5.4, for any initial data in the four-sided region which is periodic in y with period 
P, we have a uniform estimate on (u^ ,p^, p'^) as the periodic solution to the viscous 
equations (|5.2ip . From the equations in (j5.23p . we have 



V u 19^ 



dx{pu) + dyipv) = ^Ri + ^i?2 + e-^{pq) - ep{dy9f 

" " (5.32) 

B{x,y)+e-dl{pq)-ep{dye)\ 



where 

B{x, y) = ^sme( - (pq^ sin^ 6 - -)t^^^ - pq^ sm{2e)t^^^ - (pq^ cos^ 9 - -)tf-^ 
q \ p p 

+ ^ cos ( - {pq^ sin^ 9-^) t^^j - pq^ sin(20)f - (pg^ cos^ 9 - i) f 

Our L°° estimate in §5.4 guarantees that B{x,y) is uniformly bounded with respect to e. 
Using the periodicity, we have 



JO 



1 [^^ f d q 

-dl{pq)dydx = / -^dy{pq)dydx 

q ^ Jo Jo q 



Jo Jo q Jo 



dydx. 







Now integrating both sides of (j5.32p over {{x,y) : < x < xi,0 < y < P}, we find 
rxi pP ,„3{a 



Jo Jo ^ q ' 

I B{x,y)dydx- {{pu){xi,y) - {pu){xo,y))dy 
Jo 



XI fP 





<C, 

where C > is independent of e, but may depend on xi and P. This implies that 
\fedy9, \fedyq are in Lf^^{Vl) uniformly in e. 
Therefore, we have 
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Proposition 5.1. (i) Consider the viscous system (|5.2ip in 0, = {(x,y) : x > 0,y G M} 
with periodic initial data {q,6)\^=Q = (^o(y), 6*0 (?/)), then 

\/edyq, \fedy are in Lf^^{Q) uniformly in e; 

(ii) If we replace dy in system ()5.2ip by d^, the initial data {q,6)\x=Q = {qo{y),do{y)) by 

{q,O)\y=o = {qo{x),0o{x)) (5.33) 

which is periodic in x with period P and, in addition, we assume that the metric E{x) = 
G{x) is also periodic with period P, then the periodic solution with period P satisfies that 

\/edxq, \fedxQ are in L^^J^) uniformly in e, 

where Q = {{x,y) : xq < x < xi,y > 0}. 

Using Proposition 15.11 and the viscous system (j5.2ip , we conclude that 

dxM" - dyL", dxN^ - dyM"" are compact in H[^l{n). 

Since 7 G C7\ we conclude that {L%M^,N^) = ^{1" , , N") satisfies Framework (A) 
in §4. Then the compensated compactness framework (Theorem 4.1) implies that there is 
a subsequence (still labeled) {L^ , , N^){x,y) that converges weak-star to {L,M,N) as 
e ^ such that the limit (L, M, N) is a bounded, periodic weak solution to the Gauss- 
Codazzi system Therefore, {L,M,N) is a weak solution of i^M)-i^l!). We 

summarize this as Proposition 5.2. 

Proposition 5.2. For either initial value problem {i) or (ii) of Proposition \57l\ (L^,M^,N^) 
possesses a weak-star convergent subsequence which converges to a periodic weak solution of 
the associated initial value problem for the Gauss-Codazzi system (j2.6p ^ (j2.7p when e ^ 0. 

5.6. Existence of isometric immersions: Main theorem and examples. We now 

focus on the case (I5.27P : 

F = 0, and E = G depends only on x. 

to state an existence result for isometric immersions and analyze examples for this case. 
Let us look for a special solution: 

(^j?) = (0)/3) (constant state), 
for the Gauss-Codazzi system for the case (I5.27p . In this case, 

ill -2^+^' il2-U, i22 - ^11 ^12-2^ + 2^' ^22 -U, 

and the Gauss-Codazzi system (|5.ip becomes 
When q{x) = f3, this reduces to 

1 7' _ E' 
]]^~ ~ ~2E' 



or 



1 k'{x) _ E' 
~^ k{x) ~ 



(5.34) 
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Hence, q{x) = [5 becomes an exact solution precisely in this special case. Our theorem 
given below shows that, in fact, we can satisfy the prescribed initial conditions in this 
special case and that, for this choice of F, G, there exists a weak solution for arbitrary 
bounded data in our diamond-shaped region when a G (1,/?) (see Fig. 3). 

Consider the initial value problem for the Gauss-Codazzi system (j2.6p ^ (j2.7p with initial 
data 

(9, 0)U=o = (qoiy), Ooiy)), y e M, (5.35) 

or 

{q,0)\y=o = iqoix),eo{x)), x€M. (5.36) 
Our next result shows that, for this choice of E, F, G, there exists a weak solution for 
arbitrary bounded initial data in our diamond-shaped region when a G (1, /5) (see Fig. 3). 

Theorem 5.1. Assume that the initial data (j5.35p . or (j5.36p . is L°° and lies in the 
diamond- shaped region of Figs. 3-4. Then 

(i) The Gauss-Codazzi system (j2.6p - ()2.7p has a weak solution with the initial data 
{q,9)\x=o = iQoiy),do{y)). This case includes Example 5.1 for the catenoid with the metric 



E{x) = G{x) = {cosh{cx))f'^-^, F{x) = 0, c ^ 0, and (3 > I, and Example 5.2 for the 
helicoid (in isothermal coordinates) with E{x) = G{x) = ^A^ -|- ^(A^e^^ -|- e~^^), F{x) = 0, 
and (3 = \f2. 

(ii) The Gauss-Codazzi system ()2.6p - (|2.7p has a weak solution with the initial data 
{q,9)\y=o = {Qoix),OQ{x)). This case includes the catenoid with the metric E{x) = G{x) = 
{cosh{cx))'^^^ ^""^^ F{x) = 0, c 7^ 0, and /3 > 1, and the helicoid (in isothermal coordinates) 
with E{x) = G{x) = + ^A^e^^ + e'^^), F{x) = 0, and (3 = V2. 

Proof. We start with case (i). 

Step 1. For the initial data {qo{y),9o{y)), we can find {qg (y) , Oq (y)) for P > such 
that 

(i) qQ ,0q G C^(M), qQ , 9^ are periodic with period P; 

(ii) q^ — > qo, 9q 9o a.e. in M and weakly in L°°(]R) as P — > cxd. 

In particular, the functions qQ and ^q" are bounded in L°°(M), and {q^,9^) converges to 
{qo,0o) in Lfoc(M), p G [l,oo), as P ^ oo. 

This can be achieved by the standard symmetric mollification procedure: First truncate 
the initial data (PF_^0)^+,o) = (^-(^O) ^o)) W^+(90i ^o)) in the interval — < ?/ < -f 
and make the periodic extension to the whole space y G M, and then take the standard 
symmetric mollification approximation to get the G°° approximate sequence {W^q, ^+ o) 
that yields the corresponding C°° approximate sequence 

(9o '^0 ) = ((cos( ^ )) , ^ ) 

converging to {qo,9o) = ((cos( ^+'"^^-''-'> ))~^, ^+.o+'^-.o ^ M as P ^ cx). Since the 

standard symmetric mollification is an average-smoothing operator, then the approximate 
sequence [q^ ,9^) still lies in our diamond-shaped region of Figs. 3-4. 

Step 2. Following the arguments in Section 5.4, we can established the uniform L°° 
apriori estimates for the corresponding viscous solutions in two parameters e > and 
P > 0. For fixed P, then we can show that there exists a unique periodic viscous solution 
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with period P to the parabohc system (j5.2ip . which can be achieved by combining the 
standard local existence theorem with the L°° estimates. The i?;~^-compactness follows 
from the argument in ^5.51 

For fixed P, letting the viscous coefficient e tend 0, we employ Proposition 5.2 to obtain 
the global periodic weak solution {L^ , , N^) of the Gauss-Codazzi system (12.6j) - (l2.7j) . 
periodic in y with period P, in the half plane {{x,y) : x > 0, y S M}. 

Step 3. Since the sequence of periodic solutions {L^ , , N^) still stays in the invariant 
region, which yields the uniform L°° bound in P as P ^ oo. This uniform bound also 
yields the //^"^-compactness of 

d^M^ - dyL^, d^N^ - dyM^. 

Using the compensated compactness framework (Theorem 4.1) again and letting P — > 
oo, we obtain a global weak solution {L,M,N) of the Gauss-Codazzi system (j2.6p - (j2.7p 
in the half plane {{x,y) : x > 0, y G M}. 

Step 4- For (ii), it also requires the periodic approximation for the metric E{x) = G{x) 
by E^{x) = G^{x) with period P in x, besides the period approximation for the initial 
data. This can be achieved as follows: First truncate the function a{x) := in the 

interval —^<x<^ and make odd extension along the lines x = — -f, respectively. 
Then make the periodic extension from the interval — < x < to the whole space x G M 

with period P and take the standard symmetric mollification approximation to get the 

p 

C°° approximate sequence a^(x) with zero mean over each period (i.e. f a^{x)dx = 0) 

that yields the corresponding C°° approximate sequence E^{x) with period P uniquely 
determined by the differential equation: 

y'(x) = a^(x)y(x), yu=o = ^(o), 

and then take the C°° approximate sequence k^{x) with period P determined by the 
differential equation: 

K'{x) = -^—^aP{x)K{x), /CU=o = k(0) 

for the corresponding /? > 1 if needed. 

Substitute = and into the viscous Gauss-Codazzi system so that the equa- 
tions in the system is periodic in x with period P. For fixed P, the same argument as that 
for the first case yields the global period weak solution {L^ , N^) of the Gauss-Codazzi 
system (I2.6|) - (|2.7p (with E{x) replaced by E^{x), periodic in x with period P, in the half 
plane {(x,y) : x G M, y > 0}. Then, using the same argument, noting the strong conver- 
gence of E^ to E, and letting P — > oo, we again obtain a global weak solution (L, M, N) 
of the Gauss-Codazzi system (|2.6p - (|2.7p in the half plane {(x,y) : x G M, y > 0}. □ 

If we repeat the similar argument for Theorem 5.1 through the corresponding vanishing 
viscosity method on the domain Q = {(x,y) : x < 0, y G M} for problem (|5.35p and the 
domain O = {(x,y) : x G M, y < 0} for problem (j5.36p . we obtain again a weak solution. 
Together, they form a weak solution in M?. As before, the associated immersion is in C^'^. 



ISOMETRIC IMMERSIONS AND COMPENSATED COMPACTNESS 



23 



Theorem 5.2. Assume that the initial data (|5.35|) . or (j5.36p . is L°° and lies in the 
diamond- shaped region of Figs. 3-4 for the case of the catenoid or helicoid metric (as 
given in Theorem 5.1), then the initial value problem (j2.6l) -( f2^ and (15.351) . or (I5.36p . 

has a weak solution in L°°(M^). This yields a C^'^{M?) immersion of the Riemannian 
manifold into M^. 

Remark 5.1. The catenoid with circular cross-section is sketched in Fig. 5. Our theorem 
asserts the existence of a C^'^-surface for the associated metric for a class of non-circular 
cross-sections prescribed at x = 0. Similarly, the C^'^-helicoid in the (x, y, z)-coordinates 
is sketched in Fig. 6. 
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